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Abstract 

Making the assumption that high energy fermions exist in the two dimensional spin- 1/2 Heisenberg antiferromagnet we 
present predictions based on the 7r-fiux ansatz for the dynamic structure factor when the antiferromagnet is subject to a 
uniform magnetic field. The main result is the presence of gapped excitations in a momentum region near (tt, tt) with energy 
lower than that at (7r,7r). This is qualitatively different from spin wave theory predictions and may be tested by experiments 
or by quantum Monte Carlo. 



PACS numbers: 75.10Jm, 71.10Pm, 74.20Mn 

It is well known in particle physics that the nature 
of elementary excitations for one and the same model 
might change as different momenta- and energy- scales 
are probed. In particular, mesons at very high energies 
are best described as pairs of fermions (quarks) while 
at low energies they are bosons. It has been suggested 
that a similar possibility exists in the Heisenberg spin- 1/2 
antiferromagnet on a square lattice, where the low en- 
ergy excitations are Goldstone bosons while the high en- 
ergy excitations show features more resembling fermions 
(loosely called spinous in the literature). Taking this as 
a real possibility we describe here the consequences of 
adding a magnetic field to the model in order to bring 
out clearly the signatures of the underlying high energy 
fermions. The predictions made in this paper can be 
tested experimentally in materials such as copper for- 
mate tetra-deuterate (CFTD) or the organic compound 
(5CAP)2CuBr4 where the exchange constant is not too 
large. They can also be tested using numerical tools such 
as quantum Monte Carlo. 

The recently observed dispersion along the antiferro- 
magnetic zone boundary (7r/2, 7r/2)-(0, tt) (lattice spac- 
ing is set to unity throughout this Letter) in experi- 
ments on CFTD[| is an indication that the nature of 
high-energy excitations in the Heisenberg model might 
be different from the low energy ones. This disper- 
sion, which also has been confirmed using Monte Carlo 
simulations , shows a broad shallow minimum around 
(0, tt). This feature which is not captured within linear 
spin wave theory (LSW) is reminiscent of a prediction 
made by Hsu||^. Treating the Heisenberg model in an 
approximation involving massive fermions, he found a 
rather deep minimum in the dispersion at (0, tt). While 
this minimum seems much to deep to explain the exper- 
iments. Ho et al.^ recently calculated the full dynamic 
structure factor using the same fermionic picture. They 
showed that the pole at (0, tt) merges with an extensive 
high-energy continuum carrying lots of spectral weight. 
For experiments this implies that the spectral peak at 
(0, tt) would be very broad and an accurate determina- 



tion of the magnon energy would be difficult, if anything, 
the high-energy continuum would make the dispersion 
appear shallower than in Hsu's original prediction. 

Neutron scattering is necessarily an indirect probe of 
spin-1/2 fermions as it measures spin-1 excitations. So in 
order to establish the existence of fermionic excitations it 
appears that one need to be able to distinguish between 
sharp peaks(magnons) and broad continua(convolution 
of fermions). However, as will be shown here one can 
bring out the difference more clearly by applying a uni- 
form magnetic field. This is known in the context of 
the Id Heisenberg antiferromagnet where the elementary 
excitations are true spinous. There the dynamic struc- 
ture factor in a magnetic field shows low energy peaks 
at incommensurate momentum- values which change with 
the magnitude of the applied magnetic field. This effect 
which has been seen experimentally follows from the 
Bethe ansatz solution, but can also be understood as a 
splitting of the mean-field spin up- and down- spinon- 
bands||]. This splitting causes new nodes to appear at 
values determined by the magnetic field thereby causing 
low-energy peaks in the dynamic structure factor result- 
ing from inter-node scattering. 

The formulation of the 2d Heisenberg model in terms 
of fermions is based on the 7r-flux ansatz Q where first 
each spin-operator is written as a pair of fermions obey- 
ing the constraint of one particle per lattice site. Then 
the resulting quartic fcrmion coupling is decoupled with 
a mean field hopping term with imaginary parts such 
that the phase change in moving a fermion around an 
elementary plaquette is tt. This mean field theory must 
be supplemented by gauge fields which enforce the single 
occupancy constraint. An alternative approach is to use 
the mean field theory as the starting point for constucting 
a variational wavefunction. In the Gutzwiller projected 
wavefunction, the doubly occupied and unoccupied states 
are removed by hand, thereby satisfying the constraint 
exactly. While the Gutzwiller projection of the 7r-flux 
state gives quite good energy, the best trial wavefunction 
is one which combines the tt-Aux state and a spin density 



wave state. The projection of this state gives long range 
Neel order with a sublattice magnetization which is in 
excellent ageement with exact computations. Thus we 
are led to consider an effective Hamiltonian whose mean 
field solution gives the 7r-flux state with the spin density 
wave. This is accomplished by introducing a weak on-site 
repulsion V. The value of V is to be determined variation- 
ally to produce the proper magnetization and should not 
be confused with the original strong Hubbard U. Adding 
also a term coupling the fermions to the magnetic field 
the resulting effective Hamiltonian reads 




where fj^ denotes creation of a fermion at site i with 
spin projection a = {— ,+} along the 5'^-axis, n^g. is 
the corresponding density operator and 0^ is a phase on 
each link on the lattice chosen such that when summed 
clockwise around an elementary plaquette the sum equals 
TT mod(27r). As the fermions carry no charge we have cou- 
pled the magnetic field as a Zeeman term to . We must 
choose S"^ (or , not S^), as couphng to cannot be 
realized experimentally, for even in a very small field the 
antiferromagnet will orient itself such that the staggered 
direction, which here is chosen to be the z-direction, is 
perpendicular to the applied field. 

The repulsive potential will be treated self-consistently. 
That is, the potential term is first decoupled writing 

<ni„ >=^+a{-iym, (2) 

where m is the staggered magnetization. The resulting 
quadratic Hamiltonian is then diagonalized and used to 
determine the left hand side of Eq. (||), which then gives 
a relation between V/ JcS, m and h, 




where the sum is taken over the antiferromagnetic Bril- 
louin zone, N is the total number of sites and Ei,E2 
are the dispersions gotten by diagonalizing the quadratic 
Hamiltonian, 

El = J^sy/Wrh/W+J^, (4) 
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where ek — (cos^ kx -h cos^ kyY^'^ ^ and a = mF/Jcff, see 
Fig. I 

To make the connection to the parameters in the orig- 
inal Heisenberg model: the exchange constant J and the 
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FIG. 1: Fermion band dispersion in the presence of a magnetic 
field along the zone diagonal around the point (7r/2,7r/2). 
Dashed line indicates branch 1 while branch 2 is solid. 

magnetic field H, we will consider the Gutzwiller pro- 
jected ground state of the quadratic Hamiltonian above 
as a variational ansatz for the Heisenberg model ground 
state. In calculating the Gutzwiller projection we pro- 
ceed as Hsu|^ and calculate the projection exactly on a 
small cluster while approximating longer range correla- 
tions using combinatorics. While Hsu used a 2 x 2 cluster 
the results quoted here is obtained from a 3 x 3 cluster. 

Expressing the spin operators as fermions, which again 
can be expressed as linear combinations of the eigenstates 
of the quadratic Hamiltonian, the Gutzwiller projected 
energy of the Heisenberg model in a fixed magnetic field 
H is minimized with respect to the parameters h and a. 
Having obtained these values V/ JeS is calculated from 
Eq. (||). This procedure does not fix the overall scale Jeff 
which indeed was taken to be the bare value J in Hsu's 
original work. Instead we will adjust Jog to fit the spin 
wave velocity as discussed below. In Fig. |^ we show the 
values of a,/i and F/ Jeff as functions of the magnetic field 
obtained in this way. We have also plotted the total en- 
ergy and the staggered and uniform magnetization after 
projection. It should be emphasized that this procedure 
of finding the effective parameters contains no adjusta- 
bles. 



In zero magnetic field the obtained value of a is 0.45 
which is close to the value 0.5 used in Ref. Q. Fur- 
thermore the energy and staggered magnetization are 
—0.333 J and 0.33 respectively which should be compared 
to — 0.334671J and 0.3070(3) obtained from quantum 
Monte Carlo simulations]^. It should be mentioned that 
the values gotten from this Gutzwiller-projected method 
is approximate in the sense that a finite-sized cluster is 
used and no extrapolation to the infinite system have 
been attempted. From Fig. it is seen that for small 
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FIG. 2: Parameters determined from the Gutzwiller projec- 
tion on a 3 X 3 cluster. The graphs show: (a) energy, (b) a, 
(c) h, (d) V/JcB, (e) Staggered magnetization and (f) Total 
magnetization as functions of the applied magnetic field 



magnetic fields a stays roughly constant while h increases 
approximately linearly with H/J. For extremely high 
magnetic fields close to the saturation field H/J = 4, the 
minimization procedure cannot alone determine values of 
a and h as there are many values of a and h maximiz- 
ing the magnetization, which is the overwhelming con- 
tribution to the energy at high fields. We have therefore 
omitted the narrow range from H/J = 3.5 and up to the 
saturation field. 

To calculate the dynamic structure factor we calculate 
the fermion pair susceptibility and employ the RPA ap- 
proximation using the values of the parameters gotten in 
Fig. H As shown in Rcf. g] this gives the correct Gold- 
stone behavior and also gives information about the full 
spectrum. At zero temperature the dynamic structure 
factor is S{uj) = —2lm.x{(jJn = u> + iS) where using the 
RPA we have 



xM = 



Xo(w„) 



1 - VXo{uJn)' 



(5) 



where all the quantities are matrices with spin and 
momentum space indices. Xoii-^n) is a short-hand for 
Xo,af3-'ysiQ,<l',^n) which is a convolution of the fermion 
propagators G/s^ and Gsa, where Gap{k,k' ,10,/) = — < 

fl3,k',uj„fa kuj ^'^'^ ^^'^ s\ivii indices. V is diagonal 
in momentum-space and the non-zero elements in spin 

space are V++; = V ;++ = — V+_;_+ = — V_+;+_ = 

V . In the presence of a staggered magnetization it is im- 
portant to note that Xo has non-zero off-diagonal Umk- 
lapp elements which will contribute to the momentum- 
space diagonal elements of x- The calculation follows 
closely that in Ref. |]lOj, and the detailed results will be 
presented elsewhere. Here we will give the qualitative 
details as well as the most important numerical results. 
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FIG. 3: Gray-level plots of the dynamic structure factor S^^ 
along the zone edges and diagonal for parameters correspond- 
ing to, from top to bottom, H/J = 0.0, 0.5, 1.0 and 3.0. The 
solid line is the LSW result, where we have used, from top to 
bottom, Jeff/ J = 1.67, 2.43, 2.87 and 3.61. 



Numerical results for the dynamic structure factor 5^^, 
that is for spin components along the magnetic field, for 
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four different values of the magnetic field are shown in 
Fig. ^ The top panel which is for the zero field case 
shows clearly the deep minimum around (0, tt) discov- 
ered by Hsu|^ as well as the merging of the high energy 
continuum with the dominant spectral peak at this point 
calculated in Ref. The thick solid line is the spin 

wave spectrum according to LSW||l^. From the second 
panel in Fig. ^ we see that for small fields the Goldstone 
mode at (tTjTt) becomes gapped pl| in accordance with 
LSW. However, unlike LSW, the spectral function splits 
into a double peak structure near (tt, tt), and the energy 
in an incommensurate ring around (tt, tt) lies slightly be- 
low that at (tTjTt). On further increasing the magnetic 
field this energy difference between the incommensurate 
ring and the peak at (tt, tt) increases as can be seen in 
the third panel in Fig. || where H/ J — 1. 

The incommensurate low-energy peaks at finite fields 
come from the 1-1-fermion scattering which also is re- 
sponsible for the remnants of the zero field minimum at 
(0, tt). For a finite magnetic field and h < 2\/2, which 
is the case almost up to saturation branch 1 develops 
characteristic energy extrema, as indicated in Fig. 
with energies ±Jcsa/2 along the contours (sin^(r2;) -|- 
sin^(rj,))^/^ — h/2, where r is the displacement from the 
(7r/2,7r/2) and symmetric points, while branch 2 keeps 
its extremal values Jcff(a^ + h?)^/"^ /2 at (7r/2,7r/2) and 
symmetric points. The fact that the energy extremal val- 
ues for branch 1 are independent of h for h < 2^/2 has 
important consequences as it implies the presence of rel- 
atively low-energy excitations ~ Jeff a due to 1-1 fermion 
scattering over almost the whole magnetic field range. 
While the energies of these excitations do not depend on 
h, they depend on H thru Jeff (a is only weakly depen- 
dent on H). 

To explain why the low energy structure around (tt, tt) 
is lower than at (tt, tt) one can consider the coherence 
factors which appears in the calculation of Xo^- These 
coherence factors Cn, C12 and C22, associated with the 
three different types of scattering, are proportional to 
(1 + g), (1 — g) and (1 -I- g) respectively, where 



cos 



i{kx) COS(kx + Qx) + COSjky) COSjky + Qy) 

For q, the momentum of the pair-excitation equal 0, 
g = 1, while it is —1 for q = (7r,7r). Because of this, 
1-1-fermion scattering will be suppressed in the region 
very close to (tt, tt) and the dominant spectral peak at 
(tt, tt) will follow the characteristic energy of 1-2-fermion 
scattering which depends on h. Away from (tt, tt), 1-1- 
scattering, which does not depend on h, reappears giving 
rise to the low energy ring-like structure. 

We have plotted the LSW result for comparison, where 
we have determined Jeff so that the velocity of the re- 



maining Goldstone boson equals the spin-wave velocity 
in LSW. As can be seen from Fig. || LSW with Jeff de- 
termined this way gives consistently a lower value for the 
gap at (tt, tt) than is gotten using the fermion picture [p^. 
However, independently of the value of Jeff, the fermion- 
result is still very different from LSW, especially the pres- 
ence of low energy excitations (lower than at (7r,7r)) in 
an incommensurate ring around (7r,7r). 

At high magnetic fields the phase space for 1-1- and 1- 
2-scattering are well separated in energy, and so for high 
fields the energy position of the dominant spectral peak 
changes very rapidly close to {tt,tt). This is clearly seen 
for H/J = 3 in the bottom panel of Fig. ^. 

We hope that the predictions presented here can be 
checked in the near future, either using neutron scatter- 
ing experiments or numerical calculations. Observation 
of the anomalous features will provide strong support for 
the fermion formulation of the Heisenberg model. 

P.A. Lee acknowledges the support from NSF grant 
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